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Abstract 

We consider the production of baryon number from collapsing domain walls, 
and in particular examine the magnitude of CP violation which is required in such 
schemes. Taking the conventional solution to the domain wall problem in the Next- 
to-Minimal Supersymmetric Standard Model as an example, we show that the ob- 
served baryon assymmetry of the universe may have been generated, even if the 
initial explicit CP violation in the Lagrangian were so small {i.e. gravitational) that 
it could never be experimentally detected. This is possible by having the explicit 
CP violation affect the way in which the walls collapse, rather than be responsible 
for the generation of baryon number directly. Net baryon number is created at the 
domain walls by the spontaneous breaking of CP. 
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1 Introduction 



In models of baryogenesis at or close to the electroweak phase transition, one of the three 
Sakharov conditions Q (no thermal equilibrium) is fulfilled by the phase transition itself, while 
the remaining two conditions (B, C and CP violation) are provided by sphaleron mediated 
processes and some extension to the Higgs sector usually involving an explicit CP violating 
phase 0, ^. So far, the study of baryon number production involving topological defects has 
mainly addressed the first of these conditions. That is, it has been shown that the departure 
from thermal equilibrium may be provided by the collapse of cosmic strings or domain walls 
|§, ^. In this paper we shall address the second. It is generally assumed that there must be 
an explicit CP violating extension to the Higgs sector similar in magnitude to that required for 
electro-weak baryogenesis, in order to bias the production of baryons. Such CP violating terms 
may eventually be detected through their contribution to the neutron electric dipole moment 
for example |Q. Here we shall show that this is not a necessary condition for baryogenesis 
from topological defects. In fact, provided that there is spontaneous breaking of CP when the 
domain walls form Q , it is possible to generate the observed baryon asymmetry with additional 
explicit CP violating terms which are gravitationally suppressed, and which will therefore never 
be detected. We stress the difference between 'spontaneous' violation of CP which is responsible 
for the local production baryon number at the domain walls, and the 'explicit' CP violation 
which is needed in order to have a global excess of baryons. In this respect, our picture is 
similar to that proposed in Ref.Q. 

Our argument can be summarised as follows. Because the domain walls in question result 
from a breaking of CP, any particular wall is not CP invariant. Global CP invariance is provided 
by the fact that there exist different types of wall which are CP conjugates of each other. The 
mechanism which is invoked in order to remove the walls need not involve large terms in the 
Higgs potential. In fact for walls in which the higgs VEV, v, is of the order the electroweak 
scale (or larger), if the degeneracy in the minima is broken by gravitational couplings of order 
v^/Mpi, the walls will disappear well before the onset of nucleosynthesis Since individual 
walls are not CP invariant, it is possible to generate a sufficient baryon number with explicit 
CP violation of order v^/Mpi. 

We shall demonstrate this using the next-to-minimal supersymmetric standard model (NMSSM). 
It should be borne in mind however that our discussion applies to any model in which the spon- 
taneous breaking of CP produces domain walls. By choosing an example with three phases we 
are perhaps making things more difficult for ourselves, since models with more than two phases 
have their own special problems (some of which will be addressed in Ref . pOt| ) . However it is 
interesting that a case with the required properties exists already in the literature. 



The NMSSM |11| is an extension of the usual minimal supersymmetric standard model 
(MSSM) |jl2|, in which the usual two Higgs doublets Hi and H2 which are necessary to give 
masses to the up and down type quarks are supplemented by a singlet Higgs superfield A^. 
The usual n term in the lagrangian, 11H1H2, is then eliminated by invoking a symmetry 
under which every chiral superfield $ transforms as $ ^ e^'^*/^<I>. The allowed terms in the 
superpotential are then XNH1H2 — ^N^, in addition to the usual fermion mass generating 
Yukawa terms, while the Higgs part of the soft supersymmetry breaking potential is extended 
by the inclusion of two more extra trilinear soft terms Ax and in place of the MSSM term 
Bjj,HiH2 to become 

Vs"o}r = -XAx{NHiH2 + h.c.) - '^Ak{N' + h.c.) 
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+ml^\Hi\^ + mljH2\^ + m%\N\^ (1) 

where H1H2 = H^H2 — H^H^, and we shall hereafter drop the index for neutral Higgses. 

The primary motivations for the NMSSM are the elimination (or at least reparametrisation) 
of the /i problem, which is that it is not clear what could be the origin of a /i parameter in 
fj,HiH2 which is of order the electroweak scale; that it allows the evasion of the usual MSSM 
Higgs mass bounds |jT^; and the fact that this relatively minor alteration to the model gives 
an extremely rich and complex Higgs and neutralino phenomenology which can be significantly 
different from that in the MSSM [|l4|. 

When electroweak symmetry breaking occurs, the three neutral CP-even Higgs scalars ac- 
quire VEVs. Using the parameters 

H2 = P2e'^' 

N = p.e''^ (2) 

it can be shown that any true minimum of the potential does not violate CP in the sense that 
the VEVs can always be made real by an appropriate field redefinition, up to the existence of 



three degenerate vacua related to each other by Z3 transformations |15], and hence we have 
minima with 9i = 62 = Ox = ^^j^ for integer n, Q and with pi = fj, p^ = x. Note that although 
we have imposed that the fj be real, one or two of them may still be negative. We shall refer to 
the three minima 

ei + 92 = 0,27r/3,47r/3 

Ox = 0,4^3,2^3, (3) 

as A, B, C respectively, and for convenience will assume that the evolution of the universe will 
ultimately end with phase A dominating. 

After the electroweak phase transition the universe will be divided up into regions of different 
minima separated by domain walls. In each of the three degenerate minima, there is an operation 
which performs a CP transformation in the effective low energy theorem, and which maps the 
vacuum into itself. If we define the Z3 operation tohe Z-^ : A ^ B, : B ^ C, and Z3 : C — > A, 
then the transformations are 

CPa = CP 
CPb = CPZs 

CPc = CPZ|, (4) 

where here, CP is the transformation in the full theory. In each minimum, the two false vacua 
are CP conjugates of each other. Alternatively, we could have performed a field redefinition 
such that the true minimum (here A) is CP invariant. 



2 Domain Walls in the NMSSM 

Domain walls are one of the simplest types of topological defects |^, and form whenever the 
theory in question has a discrete number of degenerate vacua, usually due to the spontaneous 



^We can use weak hypercharge to select any phase for 9i — 
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breaking of a discrete symmetry. The simplest example occurs in the case of a real scalar field ^ 
with potential V = A(0^ — z/^)^. This potential clearly has two degenerate minima with cj) = ±iy 
as a result of the Z2 symmetry — > — 0. If we look for time independent solutions of the field 
equations which are translation invariant in two of the three space dimensions, and which obey 
(f){z = —00) = —v and (f){z = 00) = f we find a solution (f) = tanh(z/5z). This is a domain 
wall, whose thickness Sz is given by 5z = (\/2Az/)~^, and which has a surface energy a given by 
more complicated models, it is no longer possible to solve the field equations 
analytically, but it is straightforward to solve them numerically. We find that 6z ~ and 
a ^ as before, where u is now some typical VEV of one of the fields, and the structure is in 
general similar. 

Turning specifically to the NMSSM, the potential for the neutral scalars takes the form (at 
tree-level) 

V = M±M) (|j^^|2 _ |^2|2)2 ^ x'\N\\m' + m') 

o 

+X'^\Hi\'^\H2\^ - Xk{N^HiH2 + N^HiHi) 
+A;2|iV|4 + \H^\^ + m^^Ji/aP + m\\N\^ 

-\Ax{NH^H2 + NH,H2) - '^Ak{N^ + N^) (5) 

With real VEVs < pi >= z^i, < P2 >= 1^2, < Px >= x our inputs are then tan/3 = ^) ^ = 
A, k, A\, Ak, while u'^ = + z^l is derived from the requirement that we have the correct Z 
mass. We choose to specify the VEVs as input parameters rather than the masses appearing in 
the potential for convenience, since we may immediately calculate the soft masses rriff^, rriff^, 
m\ from the VEVs by using the minimisation conditions. Of course, this model typically has 
several different minima, usually including for example minima with only one of the three VEVs 
non-zero, and in order to study the vacuum structure it is necessary to find all of them to ensure 
that the minimum which we are analysing is indeed the deepest one. 

Let us now turn to domain wall solutions of the field equations. These reduce to the six 
equations 

where 0i is the real or imaginary part of one or other of the three scalar Higgs fields. We 
may then impose the boundary conditions that {Hi, H2, N) are {yi,V2^x) at 2; = —00 and 
(z^ie^'^*/^, z^2e^'^*^'^, a;e^'^*/^) at z = 00. It is straightforward to find solutions to such equations 
numerically, and by appropriate Z3 field redefinitions it is clear that walls with the same structure 
exist between any two pairs of vacua. 

A typical solution is shown in Figure 1, where we show the absolute values, phases, and 
energy density as a function of z in the wall region. The input parameters are \ = k = 0.2, 
Ax = Ak = lOOGeV, tan/3 = r = 2. Here the total surface energy density of the wall is 
7.1 X 10*5GeV^ after we have subtracted the vacuum energy density, while the wall thickness is 
around 0.02GeV~^, in reasonable agreement with the approximate arguments given above. It 
should be noted that even in the centre of the wall the VEVs are not zero, and so electroweak 
symmetry is not restored. 

In fact, as the parameters are varied a very wide range of different behaviours and structures 
for the wall can be seen, with typically (for tan/3 > 1 and r > 1) Px remaining large over 
much of the region, while P2 > always but may become quite small near the centre of the 
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wall. The phase behaviour shown in Figure lb, where the U{1)y phase 9i — 62 goes from 27r 
to continuously across the range is not universal but is typical. Unlike that for all the other 
variables, most of this change in 61 — O2 is outside the wall region, but we have explicitly checked 
that changing the size of the box does not have any significant impact on the total energy or 
the shape of the field configuration. 

An example of a set of parameters for which electroweak symmetry is virtually restored is 
shown in Figure 2. Here \ = k = 0.1, Ax = Ak = 250GeV, tan/3 = 2, r = 5. The total wall 
energy is 2.2 x 10''GeV'^, rather higher than before because the singlet VEV is larger, while the 
wall is now slightly wider. Although only p\ is ever zero inside the wall, p2 falls to under 3GeV, 
and is less than lOGeV for a region of width ~ 0.04GeV~^. 

Of course, these are just two of a multitude of possible sets of parameters, each of which will 
give a wall with possibly very different characteristics. We also remark that there can be more 
than one solution for a given set of parameters, although for those shown there are no other wall 
solutions with higher surface energy. Since the VEVs of the higgs fields do not go through the 
origin, there is the possibility for 'double' (and also triple in this case) wall systems with the 
same phase on either side to form in the manner described in Ref.|16|]. We shall assume that 
this does not occur, or at least that if it does the double walls are unstable to the formation of 
holes by quantum tunneling, which then expand under the surface tension destroying the wall. 
Our primary conclusion must be that for at least some sets of parameters, the domain walls 
possess exactly the properties which we will require in order to have them driving baryogenesis. 



3 Baryogenesis from Wall Networks 

Having established this fact, let us go on to examine the possibilities for baryogenesis. After the 
phase transition we have an 'emulsion' of three phases separated by highly convoluted domain 
walls. CP is also spontaneously broken by the phase transition, but as yet there is no explicit 
CP violation. In fact, as we have seen, when going through a domain wall from A ^ C, B ^ A 
or C ^ B, the phase changes of the Higgs fields are equal and opposite to those occurring 
when going from A^B,B^CoiC^A. We shall refer to these two types of transition as 
'positive' and 'negative' respectively. The walls are not invariant under CP, and inside them the 
electroweak symmetry is restored if the vacuum expectation values of the Higgs fields vanishes 
(which, as we have seen, may or may not be the case depending on the details of the Higgs 
sector). Thus we shall assume that baryon number violating transitions will be in equilibrium 
in these regions, at plasma temperatures close to the phase transition. As domain walls move 
through space, the time-dependent change of phase of the Higgs fields occurring inside the walls 
will give rise to a non-zero chemical potential for baryon number and baryons will therefore be 
created. 

Cosmology dictates that there is some mechanism which removes the walls and one sug- 
gestion, originally by Zel'dovich et al |^], is that the degeneracy of the vacua may be slightly 
broken, eventually leading to the dominance of the true vacuum. This point of view was recently 
supported by Rai and Senjanovic |17], who argue that gravitational interactions may explicitly 



break the discrete symmetries causing a slight non-degenaracy in the minima of the Higgs of or- 
der e ~ /Mpi (where u is a generic Higgs VEV of order M^y in this example). This suggestion 



was applied to this particular model in the context of string theories by Ellis et al |18]. 

We should point out two possible problems with this solution to the domain wall problem for 
this particular model. The first is the problem of destablising divergences which may generate 
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a large VEV for the singlet, and so destroy the solution which supersymmetry provides for the 
hierachy problem |19]. This is a potential fault in any model which includes gauge singlets. It 
is not clear which operators may be generated at the Planck scale or with what coefficients, 
but we note that the CP-violating gravitationally-suppressed operators which are necessary 
to remove domain walls do not in themselves generate such large singlet VEVs. Connected 
with this problem is the fact that if we break the Z3 symmetry even by gravitational terms, we 
reintroduce the /x-problem since without the Z3 symmetry there is nothing to prevent ^ becoming 
large. These points detract from the NMSSM but are unavoidable; unless we are prepared to 
complicate the model by invoking inflation with reheating to a temperature less than the weak 
scale (and probably the Affieck-Dine mechanism for baryogenesis) , we must certainly break the 
symmetry explicitly. These are problems for the NMSSM as a whole and are secondary to our 
present more general aim of showing that domain walls can induce baryogenesis with small CP 
violation. We will not discuss them further, but will simply bear in mind that a full resolution 
of these problems of the NMSSM seems to require a greater understanding of the structure at 
the Planck scale. 

The removal of the false vacua (and therefore the domain walls) proceeds as follows. For 
friction-free motion, the typical curvature scale, R, of the wall structure evolves roughly as 
the time for models with Zjsi symmetry. Since we are not interested in the precise power law 
behaviour of the curvature scale, we shall neglect the conformal stretching due to the expansion 



of the universe. For detailed discussions of these points see Refs. [2C , 21 1 . We shall also neglect the 



effects of friction on the motion of the walls. In fact this may be important at lower temperatures 
for domain walls associated with higgs fields. This is due to the walls' interaction with particles 
in the plasma, most importantly the bottom quarks, which are reflected off them with probability 
proportional to m^/p^ where p is the particle's momentum. Thus friction is unimportant for 
temperatures between Ew and 10 GeV. When the motion of the walls is friction dominated, 
they reach a terminal velocity determined by their curvature and by the density of the plasma. 
The typical curvature scale of the walls then increases as t^^'^ rather than t [^]. These points 
will be discussed in detail in [10|. 

Once the curvature scale has exceeded a critical value, i.e. when the pressure dominates 
over the tension, e > a/ R where a is the surface energy density, the domains of true vacuum 
begin to dominate and expand into the two domains of false vacuum. However, since CP is only 
broken spontaneously, any mechanism which removes the walls generates as much matter as 
anti-matter. In this case, the true vacuum. A, invades an equal area of B and C when it finally 
dominates {B and C must be degenerate if CP is not explicitly broken), and the production of 
baryons from negative walls exactly cancels that from positive ones. 

Spontaneous CP violation per se is therefore not enough to generate a net baryon number. 
What is also required is some additional explicit CP violation in the Lagrangian, and this is 
where this paper differs from previous discussions. Previously attention has almost always been 
focussed on the biasing of the baryon number production directly at the collapsing domain walls 
(an exception being the scenario examined in Ref.[Q, which bears some resemblances to this 
picture). Thus any explicit CP violation that was added to the Lagrangian was incorporated 
linearly into the production of baryon number. The resultant models required relatively large 
CP violating phases in the higgs sector. 

However, even tiny (of order v^/Mpi) CP violating terms will clearly effect the way the 
domain walls collapse, and, as argued in Ref.|17], there is no reason why gravitational terms 
that break the Z3 should not also break CP. What we propose therefore, is that no two of the 
vacua are degenerate, so that C has a higher vacuum energy than B, which has a higher vacuum 
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energy than A. As the wahs collapse therefore, A domains will invade both B and C, B domains 
will invade C but be invaded by A, and C will be invaded by both A and B. In order to show 
that this can generate a significant baryon asymmetry, consider the extreme case, in which the 
C phase has 'much' higher vacuum energy than B. Then the first pressure driven process to 
operate as the scale of the wall network increases is the collapse of C domains to be replaced 
by A and B. Both positive and negative walls will quickly accelerate to the speed of light, and 
the net baryon number generated will be close to zero. Now there will be only A and B phases 
left. The remaining B is finally removed when the non-degeneracy in A and B vacua becomes 
dominant. But the only walls which can do this are the positive B ^ A ones, and so there is 
clearly the potential for generating net baryons. 



4 Simulation of the Wall Networks 

The number which we need in order to be able to estimate the baryon number is the average 
number of positive and negative walls which pass through a given point during the whole process. 
In order to show that this number can be close to one, we have simulated a Z3 domain wall 
network evolving in 2 dimensions, in Minkowski space (i.e. neglecting the conformal stretching). 
We did this following Kawano ||2^. First we begin with an arbitrary distribution of the three 
phases. The probability for each phase is Pa = Pb = Pc = 0.33 which interestingly is barely 
enough for them to percolate in three dimensions. Simulations on a cubic lattice give the 
percolation threshold to be 0.31 and simulations in continuum percolation theories give 0.295 it 



0.02 j^. Thus the structure is expected to be tenuous and highly convoluted {i.e. 'spaghetti'- 
like). The walls are then divided into small lengths and released from the (in this case square) 
lattice. The evolution at each time step is determined by applying the equations of motion 
locally, taking the mass (proportional to the length) of the walls to be concentrated at the 
vertices between straight sections. In this we differ slightly from Kawano, who calculated the 
local curvature, since this enabled us to treat vertices with two and three walls attached on the 
same footing. In addition we did not include toroidal boundary conditions but let the ends of 
the walls slide along the edge of the box. We therefore do not expect our results to be accurate 
when the curvature scale is of the same order as the size of the box. Our basic unit for the 
simulation is shown in figure 3. A typical point, ro, is connected to up to three other points. 
Each line has a perpendicular vector e^j associated with it which describes the magnitude and 
direction of the pressure acting on it. The rest mass of the vertex is given by half the sum of the 
lengths multiplied by the surface density u; 

"^0 = ~'"o|- (7) 

i 

The force is given by —VE at the vertex 



dpQ ^ ( Ti-ro ^ 



ds Y V \ri - ro 
where 7 = (1 — Tq)^^/^ and s is proper time, so that the acceleration is given by 

d'^ro -2Ti-To , eoi7"^ A /Y^i 1 (n.\ 

^ = E 27 ^r37^ + ^h--ol (9) 
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It is easy to verify that the continuous case is recovered in the hmit as the size of the straight 
sections goes to zero. For example, consider the polygon made of N equal straight lengths, 
whose vertices are at r from the origin, and whose internal vacuum energy density is e. The 
equation of motion above leads to 

(1 -r^) £7-3 vr 
= r ^^°^iV' 

which is that of a cylinder of radius r in the limit N ^ 00 |^^. It is convenient to scale 
everything in terms of the initial curvature scale Rq, so that pi = ri/RQ, and r = t/Ro, so that 
eq.@ becomes, 

-1T2- = Z^ 27 - + \pi-po\] 2^\pj-po\- (11) 

dt i\ \pt-po\ cr / i 

The only free parameters in the simulation are therefore the two pressure variables, cbRo/^^ and 
ecRo/cr- All our results are presented with Rq normalised to 1cm. 

In figure 4 we can see how the network behaves without the effects of pressure. The scale 
of the structure evolves at roughly the speed of light (with about one wall per horizon) growing 
proportionally to the time. (This case, together with more general Zn cases has been examined 
in more detail by Press et al ||2l|.) The important point here is that (as remarked upon in 
Ref.||2^), without pressure, the evolution of the walls is mostly a question of topology. Those 
regions which are connected to two or four external walls tends to collapse, while those which are 
connected to eight or more external legs expand, due to the tension of the external legs pulling 
outwards. One can see this by considering any three leg vertex. A three leg vertex minimises 
its wall energy by trying to adopt a position in which the angles between the legs are equal and 
120°. Squares with four external legs do this by the vertices falling inwards trying to increase the 
internal 90° angle. Octagons with eight external legs expand, trying to decrease a 135° angle. 
Hexagonal structures (i.e. honeycombs) are stable. In fact for a general sided polygon with 
N external legs, the equation of motion is easily found to be 



g-J-2 j.y 




= 1-TT7— TTt)+T7^ (12) 



where we have normalised r = R/Rq, where R is the perpendicular distance from the centre to 
the edges of the polygon, and r = t/ Rq. e is the difference in vacuum energy between the inside 
and outside of the polygon, and Rq is its initial size. 

We now introduce pressure by switching on the e above. This becomes dominant over the 
tension when 

1 (13) 



which for typical values of vacuum expectation values happens for 1cm < Rq < Im in cases 
where e is induced gravitationally. The evolution of the system with pressure is shown in figure 
5, where we have taken = 0,0.25,0.5 for the phases A,B,C respectively. As mentioned 
earlier, the evolution of the network in terms of r = R/Rq is the same for constant eRq. As in 
the Z2 case, larger structures are affected much more by the pressure. Since the structure of the 
walls is always increasing, once the pressure becomes dominant, collapse happens very quickly 
(between lO"^'' and 10~^ seconds). Thus provided that the walls do not dominate in density 
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before they collapse, and also that the entropy released into the plasma is properly thermalised 
(both of which we shall assume), there is no danger of disturbing nucleosynthesis which begins 
at ~ 1 second. 

The ratio of area cleared by positive walls to the total area we find to be 0.6, so that definining 
KBG as 

area of positive transitions - area of negative transitions , 

l^BG = (14) 

total area 

the global production of baryons is kbg ~ 0.2 of what it would be for maximal CP violation. 
For values of ec much larger than this, kbg rapidly approaches 1. 



5 Discussion. 

From this point, the analysis closely follows that in Ref.[^, and of course the same caveats apply. 
That is, we assume that the wall thickness is large enough to allow anomalous processes to occur. 
These may take the form of short range interactions of typical size {g'^T)~^, where the electroweak 
symmetry is completely restored, or if the temperature is close to the phase transition one would 
expect sphaleron-like configurations straddling the domain wall to be possible. (Ideally one 
would like to be able to find these by constructing a non-contractible loop around the domain 
wall background.) Assuming that the sphaleron rate inside the walls is 

Tb ~ KiawT)"^ (15) 
the final production of baryons is given by 

^ ~ KKBGg^^atv^i^i + ^2) (16) 

where g^: is the number of relativistic degrees of freedom contributing to the entropy, and A(6i + 
62) is the phase change of the higgs fields which is orthogonal to the Goldstone mode. This factor 
is to be multiplied by a coefficent of order 1 (see Ref.[^). Note that in contrast to 'spontaneous 
baryogenesis' scenarios in which there is a possible suppression by a factor m^/T^ [^], here 
there is no suppression, because the relevant equilibrium densities to use are the ones obtaining 
outside the domain wall where mt has its physical value. Notice that the production of baryon 
number is due to the translation of the walls, which contrasts with the case of electroweak 
strings, in which it is due to a decrease in the total volume covered by strings as they collapse; 
because of this there is no volume suppression factor (SF). The change in higgs phase is 

A{ei + 62) ~ 27r/3, (17) 

so that 

!^^w-\KBGg:'- (18) 

s 

Bearing in mind our earlier discussion, there is the possibility of a much larger biasing of 
the potential. We note that this mechanism works for more general e, provided firstly that the 
amount of explicit CP violation is of the same order as the explicit Z3 breaking, and secondly 
that they are both not so large that the walls collapse immediately on forming. Thus we are in 
the novel position of being able to place (albeit extremely weak) lower and upper bounds on the 
amount of explicit CP violation allowed. The mechanism works only when the scale at which 
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the pressure dominates is larger than the size of the protodomains (~ (si^) ^) during the phase 
transition which gives, 

e ^ glT.a. (19) 

For typical values of a we find, 

e ~ lO^GeV^ (20) 

which not surprisingly is just less than My^. In addition we require that the temperature be 
close to the weak scale, for the anomalous processes to be operative inside the wall. How- 
ever the domain sizes grow at speeds comparable to the speed of light. FRW cosmology gives 
t =2.42 g~^''^{T/M(N)-'^ sees. For T ~ Mw we find t ~ 10"^° sees. Thus it is possible for 1 cm 
size structures to grow at temperatures close to the weak scale, implying that even CP violation 
induced by gravity could be the driving force behind baryon production for weak scale phase 
transitions. 

In fact the following exercise is instructive. Suppose that the anomalous processes are 
effective down to a temperature < T^, and that the pressure and surface energy terms are 
given by 

e = v^/M ■,a = v^ (21) 

where v = 0{My/) is the VEV of the higgs fields, and M is the mass scale of the physics which 
is responsible for the CP violation. Suppose also that the curvature scale increases at some 
sizeable fraction, /?, of the speed of light, R = (3{t — tc)- Then in order for this mechanism to 
work, we require that the pressure is dominant over the surface tension for t = t^, 

e > a/P{U - tc). (22) 

This gives an upper bound on M, 

M ^ M,i (^0.3/35*-'/' (^)' (l - T^/T^)^ . (23) 

So unless is extremely close to Tc, gravitational couplings could be responsible for this mech- 
anism, giving a lower bound on e of, 

e ~ 10-^GeV^ (24) 

When this bound is saturated, on dimensional grounds one expects the contribution to the 
electric dipole moment of the explicit CP violation, to be of the order of Sdn < 10~^^ecm 

Finally, we note that this mechanism is possible for any model with a spontaneous CP 
breaking transition occuring at an energy scale, v, which is higher than the electroweak scale, 
provided that some domain walls remain at the time of the electroweak transition. In this case 
anomalous processes are guaranteed to be in equilibrium when the wall network collapses. The 
same considerations apply here. That is 

^ e ^ lO-^GeV^ (25) 

In this case the lower bound is less than what would be expected to be induced by gravity, since 
we require simply that the domain walls collapse before the electroweak phase transition whilst 
anomalous processes are still in equilibrium. 
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Figure Captions 
Figure la 

Absolute magnitures of the fields as a function of position for the parameters X = k = 0.2, 
Ax = Ak = lOOGeV, tan/3 = 2, r = 2. The three lines show, from top to bottom, p^, p2, pi- 
Figure lb 

Phases of the fields as a function of position for the same parameters as Figure la. The three 
lines show 6*+ = + 62 (solid lines) , 9^ = Oi — 62 (long dashed lines) , 6x (short dashed lines) . 
Figure Ic 

Surface energy density of the wall as a function of position for the same parameters as Figure 
la. 

Figure 2a 

Absolute magnitures of the fields as a function of position for the parameters A = /c = 0.1, 
A\ = Ak = 250GeV, tan/3 = 2, r = 5. The three lines show, from top to bottom, pa, p2, pi- 
Figure 2b 

Phases of the fields as a function of position for the same parameters as Figure 2a. The three 
lines show = 9i + 62 (solid lines) , 9- = 61 — 02 (long dashed lines) , 9x (short dashed lines) . 

Figure 3 
Basic unit for wall simulation. 

Figure 4 

Wall evolution without pressure. The four time-slices shown have time 0.5 10~^''s, 1.5 10~^'^s, 
2.5 lO^-^^s, 3.75 10~-'^°s, for upper left, upper right, lower left, lower right respectively, kbg is 
less than 0.01 always. 
Figure 5 

Wall evolution without pressure. The four time-slices shown have time 0.6 10^^'^s, 1.5 10^^'^s, 
2.4 lO^-'^'^s, 3.75 10~-'^'^s, for upper left, upper right, lower left, lower right respectively, kbg is 
-0.009, -0.004, 0.023, 0.104. Final value of kbg after all walls disappeared was around 0.2. 
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